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Managementuittreksel 


Titel 

Een  vergelijking  van  drukmetingen  met  SIP  en  de  F-PVIS 

Auteur(s) 

Ing.  M.P.  van  Rooijen  en  I.  van  der  Kolk 

Datum 

juni 1997 

Opdrachtnr. 

;  A94KL482 

Rapportnr. 

:  PML1997-A18 

In  het  kader  van  de  opdracht  A94KL482,  Schokinitiatie  Explosieve  stoffen,  werk- 
pakket  1:  ‘Detonatiekarakteristieken’  zijn  diverse  meettechnieken  voor  het  meten 
van  verschillende  detonatieparameters  geevalueerd.  Een  van  de  parameters  is  de 
detonatiedruk.  Deze  kan  met  diverse  methoden  worden  gemeten.  In  dit  onderzoek 
zijn  twee  methoden  onderzocht:  de  Shock  Induced  Polarization  (SIP)  en  de  Fabry- 
Perot  Velocity  Interferometer  System  (F-PVIS).  Van  elke  methode  is  na  enkele 
aanpassingen  de  nauwkeurigheid  bepaald. 

De  SIP-meting  berust  op  het  principe  dat  sommige  kunststoffen  worden  gepolari- 
seerd  als  er  een  schokgolf  doorheen  loopt.  De  daardoor  ontstane  polarisatiespan- 
ning  kan  worden  gemeten  en  is  een  maat  voor  de  detonatiedruk  in  het  explosief. 
Met  de  F-PVIS  kan  de  oppervlaktesnelheid  van  een  snel  bewegend  spiegelend 
oppervlak  bepaald  worden.  Door  dit  spiegelend  oppervlak  tegen  een  explosief  aan 
te  brengen  kan  met  behulp  van  de  Hugoniots  van  beide  materialen  de  detonatie¬ 
druk  bepaald  worden. 

Beide  methoden  zijn  tegelijkertijd  toegepast.  Hierbij  is  uitgegaan  van  de  bestaande 
SIP-meting.  De  F-PVIS-meting  is  tijdens  het  onderzoek  verbeterd.  De  SIP-meting 
is  een  meetmethode  om  de  schokgolfsnelheid  in  een  medium  te  bepalen,  hieruit 
wordt  de  druk  in  het  medium  berekend.  Het  voordeel  van  deze  meting  is  dat  deze 
vrij  betrouwbaar  is.  Het  nadeel  is  dat  het  aflezen  van  de  doorlooptijd  van  de 
schokgolf  door  een  SIP-plaatje  vaak  moeilijk  is.  Daar  bovenop  komt  een  extra  fout 
omdat  de  schokgolfsnelheid  omgerekend  moet  worden  naar  de  deeltjessnelheid. 
Hierdoor  kan  de  fout  in  de  berekende  druk  oplopen  tot  bijna  16%. 

Beter  is  de  methode  waarmee  direct  de  deeltjessnelheid  gemeten  kan  worden,  dit 
kan  met  de  F-PVlS.  De  F-PVIS  maakt  gebruikt  van  het  Doppler-effect:  licht  dat 
gereflecteerd  wordt  op  een  snel  bewegend  voorwerp  ondergaat  een  golflengtever- 
andering.  Met  behulp  van  een  Fabry-Perot  etalon  wordt  een  interferentiepatroon 
gegenereerd  dat  afhankelijk  is  van  deze  golflengteverandering.  Dit  patroon  wordt 
met  een  streakcamera  opgenomen.  Op  de  film  van  deze  camera  is  een  veranderend 
strepenpatroon  te  zien  waarvan  de  lengte  een  maat  is  voor  de  schokgolfsnelheid  in 
een  medium.  Hieruit  kan  de  druk  in  een  medium  met  een  maximale  fout  van 
±13,5%  bepaald  worden.  De  afstanden  tussen  de  Ujtien  zijn  een  maat  voor  de 
deeltjessnelheid.  Als  hieruit  de  druk  wordt  berekend  is  een  stap  minder  nodig  dan 


TNO-rapport 


PML  1997-A18 


bij  berekening  van  de  druk  uit  de  schokgolfsnelheid,  de  fout  in  het  resultaat  is 
hierdoor  maximaal  ±12,5%.  Er  kunnen  dus  twee  schokgolfgrootheden,  te  weten 
schokgolfsnelheid  en  deeltjessnelheid,  gemeten  worden  in  een  meting.  De  fout  in 
de  resultaten  behaald  met  de  F-PVIS  wordt  voomamelijk  veroorzaakt  door  de 
kwaliteit  van  het  beeld  op  de  film. 
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Samenvatting 


Er  zijn  diverse  meetmethoden  bekend  om  drukken  in  explosieve  stoffen  te  meten. 
Een  voorbeeld  van  een  drukmeting  die  reeds  toegepast  wordt  binnen  Divisie  Mu- 
nitietechnologie  en  Explosieveiligheid  groep  Eigenschappen  Energetische  Materia- 
len  (EM)  is  de  Shock  Induced  Polarization  (SlP)-meting.  Een  nieuwere  drukmeting 
maakt  gebruik  van  het  Fabry  Perot  Velocity  Interferometer  System  (F-PVIS).  Om 
de  nauwkeurigheid  van  deze  laatste  methode  te  bepalen  wordt  deze  vergeleken  met 
de  SIP-meting. 

De  SIP-meting  is  een  meetmethode  om  de  schokgolfsnelheid  in  een  medium  te 
bepalen,  hieruit  wordt  de  druk  in  het  medium  berekend.  Het  voordeel  van  deze 
meting  is  dat  deze  vrij  betrouwbaar  is.  Het  nadeel  is  dat  het  aflezen  van  de  door- 
looptijd  van  de  schokgolf  door  een  SIP-plaatje  vaak  moeilijk  is.  Daar  bovenop 
komt  een  extra  fout  omdat  de  schokgolfsnelheid  omgerekend  moet  worden  naar  de 
deeltjessnelheid.  Hierdoor  kan  de  fout  in  de  berekende  druk  oplopen  tot  bijna  16%. 

Beter  is  de  methode  waarmee  direct  de  deeltjessnelheid  gemeten  kan  worden,  dit 
kan  met  de  F-PVIS.  De  F-PVIS  maakt  gebruik  van  het  Doppler-effect:  licht  dat 
gereflecteerd  wordt  op  een  snel  bewegend  voorwerp  ondergaat  een  golflengtever- 
andering.  Met  behulp  van  een  Fabry-Perot  etalon  wordt  een  interferentiepatroon 
gegenereerd  dat  afhankelijk  is  van  deze  golflengteverandering.  Dit  patroon  wordt 
met  een  streakcamera  opgenomen. 

Op  de  film  van  deze  camera  is  een  veranderend  strepenpatroon  te  zien  waarvan  de 
lengte  een  maat  is  voor  de  schokgolfsnelheid  in  een  medium.  Hieruit  kan  de  druk 
in  een  medium  met  een  maximale  fout  van  ±1 3,5%  bepaald  worden. 

De  afstanden  tussen  de  lijnen  zijn  een  maat  voor  de  deeltjessnelheid.  Als  hieruit  de 
druk  wordt  berekend  is  een  stap  minder  nodig  dan  bij  berekening  van  de  druk  uit 
de  schokgolfsnelheid,  de  fout  in  het  resultaat  is  hierdoor  maximaal  ±12,5%.  Er 
kunnen  dus  twee  schokgolfgrootheden,  te  weten  schokgolfsnelheid  en  deeltjes¬ 
snelheid,  gemeten  worden  in  een  meting.  De  fout  in  de  resultaten  behaald  met  de 
F-PVIS  wordt  voomamelijk  veroorzaakt  door  de  kwaliteit  van  het  beeld  op  de  film. 
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1 0  Ondertekening . 

Bijlage: 

A  Streakopname  F-PVIS 
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1  Inleiding 


In  de  loop  van  de  tijd  zijn  er  bij  TNO-PML  diverse  metingen  ontwikkeld  om 
springstoffen  te  karakteriseren.  Dit  kan  door  bijvoorbeeld  de  schokgolfsnelheid, 
deeltjessnelheid,  druk  en  dergelijke  in  een  springstof  te  meten.  Deze  parameters 
zijn  aan  elkaar  gerelateerd  wat  beschreven  wordt  in  de  schokgolftheorie. 

Er  zijn  op  dit  moment  diverse  meetmethoden  om  de  verschillende  parameters  te 
bepalen.  De  druk  kan  bijvoorbeeld  worden  bepaald  met  de  SIP-meting  waarbij 
gebruikgemaakt  wordt  van  de  dielectrische  eigenschappen  van  kunststoffen.  Met 
behulp  van  de  schokgolf  is  bet  mogelijk  sommige  kunststoffen  elektrisch  te  polari- 
seren. 

Sinds  1994  wordt  er  met  behulp  van  licht  de  deeltjessnelheid  in  springstoffen  en 
andere  materialen  bepaald.  De  hiervoor  gebruikte  opstelling  is  de  Fabry-Perot 
Velocity  Interferometer  System  (F-PVIS).  De  meetmethode  is  gebaseerd  op  het 
Dopplereffect  waaruit  volgt  dat  gereflecteerd  licht  afkomstig  van  een  snel  bewe- 
gend  voorwerp,  bijvoorbeeld  een  spiegel  op  een  detonerend  explosief,  van  golf- 
lengte  veranderd.  De  snelheid  van  dat  bewegende  voorwerp  kan  uit  een  verande- 
rend  interferentiepatroon  worden  bepaald.  Deze  opstelling  is  in  de  loop  van  de  tijd 
verder  ontwikkeld  en  getest  door  diverse  stagiairs  [4].  Hierbij  is  altijd  met  behulp 
van  een  glasvezel  met  daarop  een  grinlens  gemeten.  Deze  grinlens  wordt  gebruikt 
om  het  inkoppelrendement  van  het  licht  aan  het  uiteinde  van  de  glasvezel  te  verho- 
gen.  Het  bewegende  voorwerp  bestaat  uit  koperfolie  of  retro-reflecterend  materiaal 
op  een  explosief.  Ondanks  goede  resultaten  zijn  nog  diverse  punten  te  verbeteren 
aan  deze  meetmethode. 

Het  probleem  met  deze  oude  en  nieuwe  metingen  is  dat  de  nauwkeurigheid  en 
reproduceerbaarheid  niet  bekend  is.  De  resultaten  zijn  overeenkomstig  theoretiscbe 
waarden  maar  eventuele  verschillen  in  de  meetwaarden  van  een  en  dezelfde  meet¬ 
methode  (reprodueeerbaarheid)  en  tussen  de  methoden  onderling  (nauwkeurig¬ 
heid)  zijn  niet  bekend.  Door  nu  meerdere  metingen  tijdens  een  experiment  uit  te 
voeren  kunnen  resultaten  worden  vergeleken  en  eventuele  verschillen  mogelijk 
worden  verklaard^  Het  onderzoek  is  grotendeels  uitgevoerd  door  I.  van  der  Kolk, 
een  afstudeerder  van  de  r.k.  T.H.  Rijswijk  [7].  Dit  rapport  is  een  samenvatting  van 
zijn  bevindingen. 


I  Dit  onderzoek  is  uitgevoerd  in  het  kader  van  de  opdracht  A94/KL/482, 
Schokinitiatie  Explosieve  Stoffen,  werkpakket  1 :  Detonatiekarakteristieken. 
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2  Het  Fabry-Perot  Interferometer  System 


Om  de  snelheid  van  een  reflecterend  oppervlak  te  kunnen  bepalen  is  het  Fabry- 
Perot  Velocity  Interferometer  System  ontwikkeld.  Dit  systeem  is  in  staat  zeer  hoge 
deeltjessnelheden,  zoals  deze  optreden  bij  explosies,  door  middel  van  veranderen- 
de  interferentiepatronen  te  meten.  De  gevonden  deeltjessnelheden  kunnen  met 
behulp  van  de  schokgolftheorie  worden  omgerekend  naar  drukken.  De  huidige  F- 
PVIS  maakt  gebruik  van  een  Fabry-Perot  etalon  met  grote  spiegels,  de  afstand 
tussen  deze  spiegels  is  instelbaar. 


2.1  De  principiele  werking  van  de  F-PVIS 

Laserlicht  afkomstig  van  een  argonlaser  wordt  via  prisma’s  door  lens  A  naar  de 
fibertip  geleid,  daar  wordt  het  licht  in  de  glasvezel  gekoppeld.  Aan  het  einde  van 
de  glasvezel  wordt  het  licht  uitgekoppeld.  Het  uitgekoppelde  licht  wordt  gereflec- 
teerd  door  het  te  meten  object  en  weer  terug  in  de  vezel  gekoppeld.  Aan  de  andere 
kant  van  de  vezel  wordt  het  licht  weer  uitgekoppeld  en  dit  valt  op  lens  A.  De 
fiberpunt  staat  in  het  brandpunt  van  deze  lens  zodat  een  evenwijdige  bundel  ge- 
vormd  wordt  van  het  door  de  fibertip  uitgestraalde  licht.  De  evenwijdige  lichtbun- 
del  valt  daama  op  cilinderlens  B.  Deze  focusseert  het  licht  in  een  richting  op  de 
spleet  van  de  eerste  spiegel  (diameter  50  mm)  van  het  Fabry-Perot  etalon.  Dit 
etalon  genereert  een  interferentiepatroon  dat  door  lens  C  op  de  streakcamera  wordt 
geprojecteerd. 


GLASVEZEL  LENS  A  CILINDERLENS  (B) 


ETALON  LENSC 


INTERFERENTIE 

PATROON 


Figuur  J :  Werking  van  de  F-PVIS, 

Indien  het  meetobject  bij  het  uiteinde  van  de  glasvezel  nu  plotseling  een  bepaalde 
snelheid  krijgt  zal  het  door  dit  meetobject  gereflecteerde  licht  door  het  Doppler- 
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effect  een  golflengteverandering  ondergaan,  wat  resulteert  in  een  veranderend 
interferentiepatroon.  Uit  dit  veranderende  patroon  is  de  snelheid  van  het  reflecte- 
rend  oppervlak  te  bepalen.  Doordat  een  cylinderlens  gebruikt  wordt  zal  een  vlak  in 
het  ringenpatroon  meer  belicht  worden  dan  de  rest  van  de  ringen. 

Dit  gedeelte  van  het  interferentiepatroon  wordt  op  de  spleet  van  een  streakcamera 
geprojecteerd.  Met  behulp  van  een  snel  ronddraaiende  spiegel  worden  de  punten 
op  een  film  afgebeeld  zodat  er  lijnen  op  de  film  verschijnen  (zie  figuur  2).  Bij  een 
bekende  rotatiesnelheid  van  de  spiegel  wordt  op  deze  manier  de  afstand  tussen  de 
punten  als  functie  van  de  tijd  zichtbaar  gemaakt  [4]. 


j- 


meet  grin-  Up  Up 


Us 


Us 


Us 


Figuur  2:  Het  resultaat  van  een  verschuiving  op  de  film. 


2.2  Het  Fabry-Perot  etalon  (FP) 

Het  Fabry-Perot  etalon  is  het  hart  van  de  F-PVIS  en  wordt  gevormd  door  twee 
planparallelle  spiegels  met  een  reflectiecoefficient  van  ongeveer  95%.  Een  bundel 
laserlicht  die  op  de  eerste  spiegel  van  het  etalon  invalt  wordt  op  deze  manier  voor 
95%  gereflecteerd  en  voor  5%  doorgelaten  in  het  etalon,  dit  is  te  weinig.  Daarom 
is  in  het  etalon  een  spleet  van  1 00  p.m  breedte  in  de  reflectiecoating  van  de  eerste 
spiegel  gemaakt  zodat  meer  licht  wordt  ingekoppeld.  De  cylinderlens  zorgt  dat 
zoveel  mogelijk  licht  op  deze  spleet  wordt  gefocusseerd.  Dit  licht  valt  direct  in  op 
de  tweede  spiegel  waar  5%  licht  uit  het  etalon  wordt  gekoppeld  en  95%  licht  wordt 
gereflecteerd.  Zo  wordt  bij  elke  reflectie  95%  van  het  licht  gereflecteerd. 
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spiegelafstand 


95%  REFLECTERENDE 
SPIEGELS 


Figuur  3:  Fabry-Perot  etalon. 

Het  licht  heeft  nu  een  faseverschil  met  het  licht  dat  bij  de  eerste  keer  wordt  uitge- 
koppeld.  Dit  faseverschil  wordt  veroorzaakt  door  de  verschillende  afstanden  die 
het  licht  aflegt  (optisch  weglengteverschil).  Indien  het  onderlinge  faseverschil  een 
veelvoud  is  van  een  hele  golflengte  plus  een  halve  golflengte  zal  het  licht  uitdoven 
Bij  een  faseverschil  dat  overeenkomt  met  een  veelvoud  van  een  hele  golflengte  zal 
het  licht  constructief  interfereren. 


Voor  het  faseverschil  in  het  uitgekoppelde  licht  geldt: 
2-n-d- cos6  „ 

(p  = - Z  -K 

A, 

waarin; 

(p  :  faseverschil  [rad]; 

n:  brekingsindex  stof  tussen  spiegels  [-]; 

d:  afstand  tussen  spiegels  [m]; 

0;  hoek  tussen  normaalspiegel  en  lichtstraal  [rad]; 
A:  golflengte  van  het  licht  [m]. 


Voor  de  intensiteit  van  het  licht  als  functie  van  het  faseverschil  geldt: 


(1 


-Rf  +4-R-sin^ 


£ 

V2 


J 


waarin: 

R:  reflectiecoefficient  van  de  spiegels  [-] 


(1) 


(2) 


Hieruit  valt  af  te  leiden  dat  een  maximum  optreedt  indien  de  sinusfactor  0  is,  wat 
voorkomt  bij  een  faseverschil  gelijk  aan  m2Tt  (meZ),  oftewel  een  hele  periode.  De 
doorgelaten  hoeveelheid  licht  is  minimaal  als  de  sinusfactor  1  is.  Het  faseverschil 
is  dan  gelijk  aan  (m2+l)-7r.  Hieruit  blijkt  dat  het  interferentiepatroon  uit  concentri- 
sche  ringen  moet  bestaan.  Als  men  achter  het  etalon  een  convergerende  lens  (lens 
C)  plaatst  zal  men  het  interferentiepatroon  kunnen  afbeelden  op  een  scherm.  De 
diameter  van  deze  ringen  wordt  dan  gegeven  door  de  formule: 
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Dj 


2-/. 


(3) 


waarin: 

D  :  diameter  ring  [m]; 

fc  :  brandpuntsafstand  lens  [m]; 

j :  ringorde  [-]; 

c:  lichtsnelheid  [m-s‘*]. 


2.3  De  tijdsresolutie  van  het  etalon 


Zoals  eerder  vemeld  is  bestaat  het  FP-etalon  uit  twee  planparallelle  spiegels  waar- 
tussen  het  licht  been  en  weer  wordt  gereflecteerd  .  Doordat  het  licht  verschillende 
afstanden  aflegt  zal  een  faseverschil  ontstaan  in  het  uitgekoppelde  licht.  Door  het 
afleggen  van  deze  verschillende  wegen  krijgt  het  licht  behalve  een  onderling 
faseverschil  ook  een  onderlinge  vertraging  die  bij  veel  reflecties  behoorlijk  op  kan 
lopen  zodat  het  uitgekoppelde  licht  uitgesmeerd  wordt  over  de  tijd. 

Bij  dynamische  metingen,  waarbij  het  interferentiepatroon  dus  moet  veranderen,  is 
dit  effect  zeer  belangrijk,  het  bepaalt  in  belangrijke  mate  de  tijdsresolutie.  De  tijd 
die  nodig  is  om  90%  van  het  ingevallen  licht  uit  te  koppelen  noemt  men  de  fill- 
time  en  wordt  berekend  met  de  formule; 


t 


90% 


2-d-ln(\-R) 
c-  In/? 


(4) 


waarin: 

t9o%:  tijd  voor  het  uitkoppelen  van  90%  ingevallen  licht  [s]. 


2.4  Snelheidsberekening  uit  interferentiepatronen 

De  concentrische  ringen  waaruit  het  interferentiepatroon  bestaat  kunnen  worden 
gebruikt  bij  diverse  metingen  waarbij  licht  beinvloed  kan  worden  (bijvoorbeeld 
brekingsindexverschillen  door  luchtdrukverandering,  frequentieverandering  door 
snelheid  en  dergelijke).  De  diameterverandering  van  de  ringen  geldt  dan  in  het 
algemeen  als  resultaat  van  de  meting,  deze  diameters  moet  dan  wel  nauwkeurig 
bekend  zijn. 

Volgens  Doppler  ondervindt  licht,  dat  is  gereflecteerd  op  een  snel  bewegend 
voorwerp,  een  golflengteverandering  ten  opzichte  van  het  invallende  licht  volgens: 

AX,  2- V 


X  c 


(5) 
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waarin: 

A>i:  golflengteverandering  [m]; 

v:  snelheid  van  het  bewegende  voorwerp  [m-s"']; 

0pm.:  een  snelheid  in  de  richting  van  de  invallende  lichtstraal  geeft  een  positieve 
golflengteverandering. 


Zoals  al  eerder  is  opgemerkt  bestaat  het  interferentiepatroon,  gegenereerd  door  de 
Fabry-Perot,  uit  een  aantal  concentrische  ringen.  Volgens  formule  3  is  de  diameter 
van  deze  ringen  onder  andere  afhankelijk  van  de  golflengte.  Bij  een  veranderende 
golflengte  door  het  Doppler-effect  is  de  diameterverandering  van  de  ringen  dus  een 
maat  voor  de  snelheidsverandering.  De  formule  hiervoor  is: 


V  = 


Dt  -  D] 


I  4- 


1  j 


X-c 

4^ 


(6) 


of 


16-// 


(7) 


waarin: 

v:  snelheid  van  een  bewegend  voorwerp  [m-s'*]; 

Dj:  diameter  van  de  niet  verschoven  orde  [m]; 

Dy:  diameter  van  de  bij  Dj-horende  verschoven  orde  [m]; 

Dj+i:  diameter  van  de  eerstvolgende,  niet  verschoven,  hogere  orde  [m]; 

fc :  brandpuntsafstand  lens  C  [m]. 

Het  veranderende  interferentiepatroon  is  dus  een  maat  voor  de  snelheidsv^ranrfer- 
ing.  Aangezien  de  beginsnelheid  bekend  is  (het  reflecterende  oppervlak  begint 
vanuit  stilstand)  is  de  snelheidsverandering  ook  gelijk  de  eindsnelheid. 


2.5  De  Fringe-constante 


Indien  een  orde  zodanig  verschuift  dat  deze  precies  op  de  volgende  orde  komt  te 
liggen,  dan  is  de  verschuiving  gelijk  aan  de  Fringe-constante.  Deze  Fringe- 
constante  is  te  berekenen  met  de  formule: 


F. 


c-X 

4^ 


(8) 


waarin: 

Fc:  Fringe-constante  [m-s'']. 
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3  De  Shock  Induced  Polarisation  (SlP)-meting 


In  dit  hoofdstuk  wordt  de  SIP-methode  beschreven  en  de  basisprincipes  van  deze 
methode.  De  SIP-meting  is  een  drukmeting  en  wordt  al  geruime  tijd  toegepast. 
Doel  hiervan  is  de  schokgolfsnelheid  te  meten  en  deze  met  behulp  van  de  schok- 
golftheorie  om  te  rekenen  naar  de  druk. 


3.1  De  principiele  werking  van  de  SIP-meting 

De  SIP-methode  berust  op  bet  principe  dat  sommige  kunststoffen  elektrisch  gepo- 
lariseerd  worden  als  er  een  schokgolf  doorheen  loopt.  De  daardoor  ontstane  polari- 
satiespanning  kan  worden  gemeten  met  een  snelle  oscilloscoop.  Dit  wordt  gedaan 
door  een  kunststofplaatje  tussen  een  aard-  en  een  meetelektrode  te  plaatsen  met 
daar  overheen  een  afsluitweerstand  (figuur  4).  Het  kunststofplaatje  bestaat  meestal 
uit  een  polymethylmethaacrylaat  (PMMA)  en  wordt  het  SIP-plaatje  genoemd.  Om 
het  elektrische  veld  dat  tijdens  de  meting  optreedt  niet  te  verstoren  wordt  er  ook 
nog  een  guardelektrode  toegevoegd.  Voordat  een  meting  met  het  plaatje  wordt 
uitgevoerd  is  het  belangrijk  de  dikte  van  het  plaatje  zo  nauwkeurig  mogelijk  op  te 
meten,  omdat  met  behulp  van  de  dikte  van  het  plaatje  de  schokgolfsnelheid  wordt 
bepaald. 


springstof  SIP-plaatje 


-guardelektrode 


-meetelektrode 


- > —  meetsignaal 

50  Q  afsluitweerstand 


aardelektrode 


Figuur  4:  de  SIP-meetopstelling. 

Het  meetsignaal  op  de  oscilloscoop  heeft  de  vorm  zoals  getoond  in  figuur  5.  Op 
het  tijdstip  to  treedt  de  schokgolf  bij  het  SIP-plaatje  in.  Op  tijdstip  tj  heeft  de 
schokgolf  de  meetelektrode  bereikt.  De  stijgtijd  van  het  signaal  ts  -  to  wordt  be¬ 
paald  door: 

•  de  vlakheid  van  de  schokgolf; 

•  de  grootte  van  de  elektroden; 

•  de  weerstand  op  het  contactvlak  aansluitdraden(SIP-pennen)/SIP-plaatje. 
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De  stijgtijd  is  dus  eigenlijk  een  maat  voor  de  kwaliteit  van  de  meting  want  hoe 
kleiner  deze  stijgtijd  is,  des  te  beter  het  tijdstip  to  bepaald  kan  worden.  Tussen  ts 
en  ti  is  een  kromme  te  zien.  Deze  kromme  wordt  veroorzaakt  door  de  condensa- 
torwerking  van  de  opstelling.  De  meetelektrode  en  de  schokgolf  vormen  de  platen 
van  de  condensator.  Het  tussenliggende  PMMA  doet  dienst  als  dielektricum.  Door 
de  polarisatie  van  het  dielektricum  gedraagt  het  geheel  zich  als  een  opladende 
condensator.  Na  het  invallen  van  de  schokgolf  bewegen  de  platen  (en  dus  niet  de 
elektroden)  met  de  schokgolfsnelheid  naar  elkaar  toe. 

Er  wordt  echter  ook  steeds  meer  PMMA  gepolariseerd  zodat  de  spanning  tussen  de 
elektroden  snel  toeneemt.  Als  de  schokgolf  de  meetelektrode  bereikt  vormen  de 
condensatorplaten  een  kortsluiting  en  gaat  de  spanning  naar  0  V. 


tijdt 


Figuur  5:  Een  voorbeeld  van  een  SIP-meetsignaal. 


3.2  Snelheidsbepaling  uit  elektrische  signalen 

De  doorlooptijd  van  de  schokgolf  door  het  SIP-plaatje  kan  worden  bepaald  uit  het 
verschil  tussen  de  tijden  to  en  t|  en  de  dikte  van  het  SIP-plaatje 


At 

waarin: 

Lsip:  dikte  van  het  SIP-plaatje  [m]; 

At;  doorlooptijd  van  de  schokgolf  door  het  SIP-plaatje  [t]. 

De  druk  in  de  springstof  kan  nu  met  behulp  van  de  schokgolftheorie  worden  bere- 
kend. 
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4  Uitvoering  van  gecombineerde  drukmetingen 


Als  men  de  Hugoniot-vergelijking  van  een  medium  wil  bepalen  moet  er  een  serie 
metingen  worden  uitgevoerd  waarin  telkens  twee  grootheden  uit  de  schokgolfthe- 
orie  tegelijkertijd  worden  bepaald.  Bijvoorbeeld  de  deeltjessnelheid  met  de  F- 
PVIS  en  de  schokgolfsnelheid  volgens  de  SIP-methode.  Met  behulp  van  de  compu¬ 
ter  kan  dan  de  relatie  tussen  de  gemeten  grootheden  worden  gezocht  (‘curve- 
fitten’).  Uit  de  fit  zijn  dan  de  parameters  voor  de  Hugoniot-vergelijking  voor  dat 
materiaal  af  te  leiden.  De  nauwkeurigheid  en  betrouwbaarheid  van  de  meetmetho- 
den  moet  dan  wel  precies  bekend  zijn. 

Als  de  Hugoniot-vergelijking  van  een  materiaal  al  bekend  is  kan  met  het  meten  van 
een  van  de  schokgolfgrootheden  worden  volstaan,  de  andere  grootheden  kunnen 
dan  berekend  worden.  Daarom  gebmiken  we,  om  de  verschillende  methoden  te 
kunnen  testen  en  vergelijken,  alleen  materialen  waarvan  de  Hugoniot-vergelijkin- 
gen  bekend  zijn. 

Om  de  diverse  meetmethoden  goed  te  kunnen  vergelijken  moeten  ze  tegelijkertijd 
in  een  meting  worden  uitgevoerd;  eventuele  spreiding  in  de  eigenscbappen  van  de 
springstof  worden  zo  geelimineerd.  Aangezien  het  oppervlak  van  het  explosief 
waaraan  gemeten  wordt  beperkt  is,  moet  er  eerst  een  goede  opstelling  worden 
gemaakt  die  eenvoudig  en  compact  is.  Omdat  ook  de  reproduceerbaarheid  wordt 
gemeten  is  het  van  belang  de  inleiding  ook  zo  reproduceerbaar  mogelijk  te  laten 
plaatsvinden.  Hierbij  wordt  gebruikgemaakt  van  een  Plane  Wave  Lens  (PWL) 
donorlading  waarop  een  TNT  booster  is  aangebracht  (diameter  60  mm;  hoogte 
20  mm).  De  specificaties  van  de  PWL  zijn: 

diameter:  60  mm; 

formule  lensbolling:  d  =  18,00  -  0,0189r2; 

d:  hoogte  lens  [mm]; 
r:  straal  lens  [mm]; 
explosief:  PBX:  RU-81. 

De  initiatie  vindt  plaats  met  een  PL  438  slagpijpje.  Tussen  het  slagpijpje  en  de 
PBX  van  de  PWL  wordt  kneedspringstof  geplaatst  (DM12)  om  een  goed  contact 
tussen  beide  onderdelen  te  verkrijgen.  Het  hele  donorsysteem  wordt  getoond  in 
figuur  6. 
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TNT 

plane  wave  lens 


PBX 


kneedspringstof 


PVC  slagpijphouder 
slagpijpje 


Figuur  6:  De  PWL-donorlading. 

Het  meetgedeelte  is  gedurende  het  onderzoek  aangepast,  met  de  voorwaarde  dat 
beide  meetmethoden  (F-PVIS  en  SIP)  tegelijkertijd  kunnen  worden  gebruikt.  De 
verschillende  configuraties  worden  in  de  volgende  paragraien  besproken,  waama 
telkens  ook  de  resultaten  met  deze  configuratie  worden  besproken. 

De  resultaten  worden  verkregen  door  meetgegevens  van  de  F-PVIS-  en  SIP-meting 
om  te  rekenen  naar  de  druk  in  de  springstof.  In  de  hiervoor  gebruikte  Hugoniot- 
vergelijkingen,  waarbij  de  parameters  met  behulp  van  een  ‘curve-fit’  bepaald  zijn, 
moeten  de  eenheden  in  km  s  '  en  GPa  worden  uitgedrukt.  Om  verwarring  te  voor- 
komen  zijn  in  dit  hoofdstuk  alle  grootheden  volgens  deze  eenheden  ingevuld.  In 
het  hieropvolgende  hoofdstuk  is  een  foutberekening  uitgevoerd  over  de  gevonden 
resultaten.  In  hoofdstuk  6  wordt  de  drukberekening  nog  eens  herhaald  maar  nu 
wordt  er  van  uitgegaan  dat  de  Hugoniot  van  de  springstof  onbekend  is.  In  hoofd¬ 
stuk  7  wordt  een  nieuwe  methode,  de  glasfiberprobe,  besproken. 
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4.1  De  originele  opstelling  (F-PSIPl) 


zilververf 


Figuur  7:  De  opstelling  bij  meting  F-PVISl. 


4.1.1  De  SIP-meting  (schokgolfsnelheidsmeting) 

Om  eerdere  meetresultaten  van  de  SIP-meting  te  kunnen  gebruiken,  is  in  eerste 
instantie  geprobeerd  zoveel  mogelijk  de  originele  meetopstelling  te  behouden. 

Bij  de  originele  SIP-meting  wordt  de  gehele  PWL-donorlading  in  aluminiumfolie 
gewikkeld  (huishoudfolie).  Op  het  TNT  wordt  een  dunner  stuk  aluminiumfolie 
(10pm)  gelegd.  Dit  folie  is  van  tevoren  gladgewreven  zodat  het  goed  contact  maakt 
met  de  springstof.  Dit  aluminiumfolie  dient  als  aardelektrode  voor  de  SIP-meting. 
Op  dit  aluminiumfolie  komt  het  SIP-plaatje  te  liggen  (zie  figuur  7).  De  dikte  van 
dit  plaatje  moet  zo  precies  mogelijk  worden  opgemeten  omdat  met  behulp  van 
deze  dikte  schokgolfsnelheidsberekening  plaatsvindt. 

De  guard-  en  de  meetelektrode  worden  gemaakt  door  de  bovenzijde  van  het  SIP- 
plaatje  met  zilververf  te  beschilderen  (dit  geleidt  de  elektrische  stroom  goed).  In 
het  midden  van  het  plaatje  wordt  in  de  zilververflaag  met  behulp  van  een  kraspas- 
ser  een  cirkel  gekrast.  De  binnencirkel  (10  mm  diameter)  wordt  de  meetelektrode 
en  het  gebied  daarbuiten  de  guardelektrode. 

Voor  het  aansluiten  van  deze  elektroden  solderen  we  koperen  staafjes  aan  de 
meetdraden.  Deze  worden  in  een  plexiglas  balkje,  waarin  twee  gaatjes  zijn  ge- 
boord,  geschoven  en  het  geheel  komt  op  het  SIP-plaatje  te  staan. 

Tenslotte  verbinden  we  de  guardelektrode  en  de  aardelektrode  met  een  afsluit- 
weerstand  (50  Q).  Met  een  coaxkabel  worden  de  meetdraden  naar  de  Digital 
Signal  Analyser  (DSA)  geleid.  Met  dit  apparaat  kan  het  zeer  korte  meetsignaal 
opgenomen  worden. 
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De  resultaten  van  de  SIP-metingen  staan 

in  onderstaande  tabel  vermeld: 

Tabel  1 : 

Resultaten  SIP-meting. 

Lsip  [m] 

At[s] 

977-10-3 

155-10-9 

De  schokgolfsnelheid  Ug  kan  uit  deze  gemeten  waarde  berekend  worden.  De 
deeltjessnelheid  Up  in  het  PMMA  kan  berekend  worden  door  deze  de  gemeten 
schokgolfsnelheid  in  de  Hugoniot-vergelijking  van  PMMA  in  te  vullen,  die  weer- 
gegeven  wordt  door; 

r,  PMMA  _  ^PMMA  oPMMA  PMMA  /w 

^  S  ~  '^0  ^  ^ 

SPMMA.  1,595  [-]; 

CqPMMA.  2,500  [km-s->]; 

PqPMM A  .  1,186  [g-cm‘2] . 


Deze  berekende  deeltjessnelheid  is  gelijk  aan  de  deeltjessnelheid  op  het  grensvlak 
TNT-PMMA.  Hieruit  kan  ook  druk  in  het  PMMA  worden  berekend. 


,PMMA 


~  i^PMMA  ■  P O.PMMA  )  ■  ^ 


p.PMMA 


+  (C, 


O.PMMA 


•Po, 


PMMA  ■ 


U 


p.PMMA 


(11) 


De  Hugoniot  voor  TNT  wordt  hiervoor  zodanig  gespiegeld,  dat  de  drukken  bij  de 
deeltjessnelheid  op  het  grensvlak  gelijk  zijn.  Het  punt  waar  de  gespiegel- 

de  Hugoniot  van  TNT  snijdt  met  de  originele  Hugoniot  geeft  dan  de  deeltjessnel¬ 
heid  Up™^  in  het  TNT.  Hieruit  wordt  vervolgens  de  druk  in  het  TNT,  P^^^  bere¬ 
kend.  De  waarden  voor  de  gebruikte  Hugoniot  van  TNT  zijn: 


S™T:  2,500  [-]; 

Co™T:  2,390  [km-s->]; 

Po™T:  1,578  [g-cm-3]. 

De  gevonden  resultaten  staan  in  de  onderstaande  tabel,  de  nauwkeurigheid  van  de 
resultaten  wordt  behandeld  in  hoofdstuk  5. 


Tabel  2:  Berekening  van  de  druk  in  TNT  uit  de  SIP-meting. 


meting 

Us 

Up 

U  TNT 

Up 

pTNT 

[-1 

[km-s""'] 

[km-s"'] 

[km-s'^l 

[GPa] 

F-PSIP1 

6,300 

2,382 

2,041 

24,129 

4.1.2  De  F-PVIS-meting  (deeltjessnelheidsmeting) 

De  bedoeling  is  deeltjessnelheid  gelijktijdig  met  de  schokgolfsnelheid  te  meten 
zonder  de  originele  SIP-opstelling  ingrijpend  te  veranderen.  Een  probleem  dat 
hierbij  optreedt  is  dat  een  spiegelend  meetobject  nodig  is  dat  een  maat  is  voor  de 
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deeltjessnelheid  in  het  TNT.  Ook  moet  een  glasfiber  met  daarop  een  grinlens 
uitgelijnd  kunnen  worden  ten  opzichte  van  dit  spiegelende  meetobject. 

Als  spiegelend  meetobject  dient  het  aluminiumfolie  tussen  het  PMMA  en  het 
TNT.  De  Hugoniot-vergelijkingen  van  PMMA  en  TNT  zijn  bekend  en  hieruit  kan 
met  behulp  van  de  schokgolftheorie  de  dmk  in  het  TNT  worden  berekend.  Dat 
tussen  de  grinlens  en  het  aluminiumfolie  een  laag  PMMA  zit  is  geen  probleem 
omdat  uit  andere  onderzoeken  is  gebleken  dat  bij  een  schokgolfdoorgang  in 
PMMA  alleen  de  brekingsindex  verandert  en  geen  verduistering  optreedt  zoals  bij 
de  meeste  andere  materialen. 

Blijft  het  uitlijnprobleem:  bij  eerdere  metingen  met  de  F-PVIS  is  de  grinlens  in  een 
driepootconstructie  gezet.  Hiermee  wordt  de  meest  ideale  uitlijning  ten  opzichte 
van  een  spiegelend  meetobject  ingesteld.  In  deze  opstelling  is  daar  te  weinig  plaats 
voor.  Een  extra  gat  in  de  houder  blijkt  de  driepootconstructie  overbodig  te  maken. 
De  grinlens  wordt  eenvoudig  in  het  gat  geschoven  en  met  de  hand  gepositioneerd, 
een  dmppel  lijm  zet  het  geheel  vast. 

De  film  (zie  bijlage  A)  laat  een  duidelijk  strepenpatroon  zien  dat  aangeeft  dat  er 
genoeg  licht  door  het  aluminiumfolie  wordt  gereflecteerd.  Verder  is  te  zien  dat  de 
strepen  plotseling  lichter  worden,  op  dit  moment  heeft  de  schokgolf  het  alumini¬ 
umfolie  bereikt.  Er  zou  hier  een  verschoven  strepenpatroon  te  zien  moeten  zijn,  dit 
is  echter  niet  het  geval.  Het  lijkt  erop  dat  het  aluminiumfolie  tijdens  de  schokgolf¬ 
doorgang  zijn  reflectiviteit  heeft  verloren.  Er  kunnen  verschillende  oorzaken  zijn: 
de  lucht  tussen  het  TNT-patroon  en  het  aluminiumfolie  ioniseert  waardoor  het 
folie  verbrandt  of  het  aluminiumfolie  verliest  zijn  reflectiviteit  door  de  schokgolf¬ 
doorgang. 

Verder  is  er  ook  een  donkere  streep  dwars  over  de  film  te  zien.  Dit  wordt  veroor- 
zaakt  door  de  schokgolf  die  de  lucht  vlak  voor  de  grinlens  ioniseert.  Uit  deze  film 
zijn  geen  deeltjessnelheden  te  bepalen. 
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4.2  Reflectie  met  een  goudlaag  (F-PSIP2) 


Figuur  8:  De  opstelling  van  meting  F-PSIP2. 

Aan  de  opstelling  voor  de  SIP-meting  hoeft  niets  te  worden  veranderd,  deze  geeft 
al  een  goed  signaal.  De  opstelling  voor  de  F-PVIS-meting  moet  wel  aangepast 
worden.  Dit  kan  echter  niet  zonder  de  opstelling  voor  de  SIP-meting  bei'nvloeden. 
De  bedoeling  is  om  aan  de  onderkant  van  het  SIP-plaatje  een  goudlaag  op  de 
dampen,  die  als  spiegelend  meetobject  dient.  De  dikte  van  deze  goudlaag  ligt  in  de 
orde  grootte  van  enkele  tientallen  nanometers.  De  invloed  hiervan  op  het  resultaat 
van  de  SIP-meting  zal  vrijwel  nihil  zijn.  Verder  wordt  er  tussen  het  TNT  en  alu- 
miniumfolie  een  laagje  siliconenvet  aangebracht.  Er  zit  dan  veel  minder  lucht 
onder  het  folie  zodat  er  minder  ionisatie  van  de  lucht  optreedt.  Siliconenolie  is 
beter  maar  trekt  weg  in  de  TNT-patroon.  Tussen  het  aluminiumfolie  en  de  goud¬ 
laag  wordt  wel  siliconenolie  gebruikt.  Verder  wordt  een  dikker  SIP-plaatje  ge- 
bruikt  waardoor  er  langer  met  de  F-PVIS  gemeten  kan  worden.  Door  de  donkere 
Streep  over  de  film  van  de  vorige  meting  is  het  idee  ontstaan  dat  behalve  de  deel- 
tjessnelheid  ook  de  schokgolfsnelheid  met  de  F-PVIS  te  meten  is.  Een  verschoven 
lijnenpatroon  geeft  immers  aan  dat  de  schokgolf  het  spiegelend  meetobject  heeft 
bereikt  en  de  verschuiving  houdt  op  als  de  schokgolf  de  grinlens  op  de  glasvezel 
raakt.  Als  de  grinlens  direct  op  het  SIP-plaatje  wordt  gezet  kan  de  doorlooptijd  uit 
de  film  worden  berekend.  De  dikte  van  het  SIP-plaatje  is  bekend  en  uit  deze  twee 
gegevens  kan  eenvoudig  de  schokgolfsnelheid  berekend  worden. 
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4.2.1  De  SIP-meting  (schokgolfsnelheidsmeting) 

De  resultaten  van  de  SIP-meting  zijn; 

Tabel  3: 

Resultaten  SIP-meting. 

Lsip  [mm] 

At  [s] 

3,940 

692-1 0-9 

De  resultaten  van  de  SIP-meting  worden  op  dezelfde  manier  verwerkt  als  bij  de 
vorige  SIP-meting.  De  gevonden  resultaten  zijn  als  volgt: 

Tabel  4:  Berekening  van  de  druk  in  TNT  uit  de  schokgolfsnelheid  van  de  SIP-meting. 


meting 

Us 

Up 

U  ■'■NT 

Up 

pTNT 

[-] 

[km-s"'] 

[km-s"''] 

[km -S''] 

[GPa] 

F-PSIP2 

5,700 

2,006 

1,722 

18,184 

4.2.2  De  F-PVIS-meting  (deeltjessnelheidsmeting) 

Het  strepenpatroon  is  zeer  duidelijk;  goud  als  reflecterend  materiaal  werkt  dus 
goed.  Op  het  punt  dat  het  strepenpatroon  lichter  wordt  is  er  kort  een  verschuivend 
patroon  te  zien.  Hieruit  kan  een  deeltjessnelheid  worden  afgeleid  ,  met  een  micro¬ 
scoop  worden  de  afstanden  tussen  de  lijnen  bepaald.  Het  verschoven  strepenpa¬ 
troon  is  onder  de  microscoop  echter  maar  zeer  kort  zichtbaar.  Door  de  film  met  een 
scanner  in  te  lezen  kan  met  behulp  van  beeldbewerkingssoftware  de  film  lets 
duidelijker  worden  gemaakt,  dit  beeld  staat  in  bijlage  A. 

Hierop  is  te  zien  dat  het  verschoven  patroon  keurig  doorloopt  tot  het  einde  van  de 
meting.  Helaas  is  dit  niet  te  zien  onder  de  microscoop  zodat  het  verloop  van  de 
deeltjessnelheid  niet  in  goede  resultaten  is  uit  te  drukken.  Wei  is  te  zien  dat  de 
deeltjessnelheid  afloopt  en  waarschijnlijk  nog  snel  ook.  Waarom  de  reflectiviteit 
van  het  goud  minder  wordt  is  niet  bekend,  mogelijk  is  de  schokgolf  te  krachtig  of 
de  drukverandering  te  groot.  Waarom  de  deeltjessnelheid  zo  snel  afloopt  valt  ook 
nog  niet  te  verklaren.  Er  kan  wel  een  doorlooptijd  van  de  schokgolf  door  het  SIP- 
plaatje  uit  de  film  worden  gehaald,  die  een  maat  is  voor  de  schokgolfsnelheid, 
zodat  nu  in  een  enkele  meting  met  de  F-PVIS  twee  grootheden  bepaald  kunnen 
worden.  Deze  meting  is  nog  eens  herhaald  onder  de  naam  F-PSIP3.  Bij  deze  me¬ 
ting  is  gebruikgemaakt  van  een  elektronische  streakcamera.  Metingen  met  behulp 
van  deze  camera  zijn  nog  in  de  experimentele  fase  zodat  resultaten  hiervan  nog 
niet  beschikbaar  zijn. 
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De  lengte  van  de  meting  en  afstanden  tussen  de  lijnen  van  het  niet-verschoven  en 
het  verschoven  interferentiepatroon  op  de  film  zijn  als  volgt. 


Tabel  5:  Meetresultaten  van  de  F-PVIS  in  meting  F-PSIP2. 


orde 

D, 

[m] 

Dv 

[m] 

If 

1 

2,135-10-3 

3,180-10-3 

2 

3,345-10-3 

4,105-10-3 

3 

4,215-10-3 

4,840-10-3 

4 

4,935-10-3 

5,480-10-3 

5 

5,565-10-3 

6,065-10-3 

5,633-10-3 

Met  de  gemeten  resultaten  van  de  F-PVIS  kan  de  deeltjessnelheid  Up  met  behulp 
van  formule  7  worden  berekend.  In  tabel  6  staat  het  gemiddelde  van  de  over  de 
verschillende  ordes  gemeten  deeltjessnelheden.  Met  deze  deeltjessnelheid  kan  de 
deeltjessnelheid  in  het  TNT  direct  worden  berekend  .  Dit  is  een  stap  minder 
dan  dezelfde  berekening  met  de  schokgolfsnelheid.  Met  de  Hugoniot  van  TNT  kan 
de  druk  in  de  springstof  worden  berekend. 


Tabel  6: 

Berekening  van  de  druk  in 
meting. 

TNT  met  de  deeltjessnelheid  uit  de  F-PVIS- 

meting 

Up 

u 

Up 

pTNT 

I-] 

[km-s-^] 

[km-s-'] 

[GPa] 

F-PSIP2 

2,640 

2,259 

22,660 

We  weten  dat  de  film  met  125  ns-mm''  wordt  beschreven.  Bij  een  bekende  lengte 
van  het  verschoven  lijnenpatroon  kan  de  doorlooptijd  van  de  schokgolf  in  het  SIP- 
plaatje  worden  berekend,  en  daarmee  wordt  de  schokgolfsnelheid  bepaald.  De  druk 
in  de  springstof  uit  de  schokgolfsnelheid  wordt  op  dezelfde  manier  berekend  als 
bij  de  SIP-meting. 

Tabel  7:  Berekening  van  de  druk  in  TNT  met  de  schokgolfsnelheid  uit  de  F-PVIS- 

meting. 


meting 

Us'''"' 

Up 

U  TNT 

pTNT 

[-1 

[km-s-'l 

[km-s-'] 

[km-s"''] 

[km-s"'] 

F-PSIP2 

5,596 

1,941 

1,666 

17,235 
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4.3  Schokgolfdemping  met  PMMA  (F-PVIS4) 


zilververf 

goudlaag 

aluminiumfolie 

PMMA 


SIP-pennen 


glasfiber 

* 


■  houder 

■  PMMA 

■  siliconenolie 

■  siliconenolie 

siliconenvet 


Figuur  9:  De  opstelling  van  meting  F-PVIS4. 

De  schokgolf  in  de  vorige  meting  is  waarschijnlijk  te  krachtig  of  de  plotselinge 
drukverandering  is  te  groot  voor  de  goudlaag  zodat  deze  minder  reflecterend 
wordt.  Om  dit  te  voorkomen  wordt  een  dun  PMMA-schijfje  (1  mm)  tussen  het 
aluminiumfolie  en  het  TNT  geplaatst.  Tussen  het  PMMA-schijfje  en  het  TNT- 
patroon  wordt  een  laagje  siliconenvet  aangebracht.  Tussen  het  schijfje  en  het 
aluminiumfolie  zit  siliconenolie.  Volgens  de  schokgolftheorie  moet  de  druk  nu 
stapsgewijs  naar  zijn  eindwaarde  klimmen.  Bij  andere  experimenten  is  gebleken 
dat  de  demping  van  PMMA  ongeveer  3%-mm-'  is,  hiermee  moet  rekening  worden 
gehouden  in  het  resultaat.  De  rest  van  de  opstelling  blijft  ongewijzigd. 

4,3.1  De  SIP-meting  (schokgolfsnelheidsmeting) 

De  resultaten  van  de  SIP-meting  zijn: 

Tabel  8:  Resultaten  SIP-meting. 


LsiP 

At[s] 

0,978-10-3 

162-10-9 

Tabel  9:  Berekening  van  de  druk  in  TNT  uit  de  schokgolfsnelheid  van  de  SIP-meting. 


nneting 

Us 

Up 

u 

pTNT 

[-] 

[knn-s-''] 

[km-s-’'] 

[km-s-^l 

[GPa] 

F-PS1P4 

6,040 

2,219 

1,903 

21 ,455 
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4.3,2  De  F-PVIS-meting  (deeltjessnelheidsmeting) 

De  film  is  zeer  duidelijk.  Er  is  een  goed  verschoven  lijnenpatroon  te  zien  dat  de 
gehele  meting  zichtbaar  blijft  zodat  er  weer  een  doorlooptijd  te  bepalen  is.  De 
verschoven  lijnen  zijn  praktisch  recht  waaruit  blijkt  dat  de  schokgolf  over  1  mm 
PMMA  bij  benadering  ideaal  en  dus  ongedempt  verloopt.  Ook  de  SIP-meting  is 
goed  verlopen. 

Deze  meting  mag  dus,  als  vergelijking  tussen  SIP  en  F-PVIS,  als  geslaagd  worden 
beschouwd.  Of  het  resultaat  ondanks  de  relatief  grote  meetafstand  een  goede 
dmkmeting  is,  blijft  een  vraag. 

De  verschoven  lijnen  zijn  tot  het  einde  van  het  experiment  te  zien.  De  afstanden 
tussen  de  lijnen  zijn  aan  het  begin  en  het  einde  van  de  verschuiving  gemeten, 
hierin  zit  geen  significant  verschil.  De  afstanden  zijn  daarom  op  een  punt  bepaald. 


Tabel  10:  Meetresultaten  van  de  F-PVIS  in  meting  F-PSIP4. 


orde 

Di 

[m] 

Dv 

[m] 

kfiim 

[m] 

1 

2,240-10-3 

3,100-10-3 

2 

3,405-10-3 

4,040-10-3 

3 

4,265-10-3 

4,775-10-3 

4 

4,980-10-3 

5,420-10-3 

5 

5,595-10-3 

5,995-10-3 

1,290-10-3 

Tabel  11: 

Berekening  van 
meting. 

de  druk  in  TNT  uit  de  deeltjessnelheid 

van  de  F-PVIS- 

meting 

Up 

u 

pTNT 

[-1 

[km-s-^] 

[km-s-’l 

[GPa] 

F-PSIP4 

2,154 

1,847 

20,424 
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label  12: 

Berekening 

meting. 

van  de  druk  in  TNT  uit  de  schokgolfsnelheid  van 

de  F-PVIS- 

meting 

Us'''"" 

Up 

U  ■'■NT 

Up 

pTNT 

[-] 

[km-s'^ 

[km-s‘^] 

[km-s"'] 

[km-s"'] 

F-PSIP4 

6,065 

2,235 

1,916 

21,705 

4.4  Mislukte  meting  dichter  op  de  springstof  (F-PSIP5) 


zilververf 

goudlaag 

molybdeenfolie 

aluminiumfolie 


Figmr  W:  De  opstelling  van  meting  F-PSIP5. 

Met  de  opstelling  van  meting  F-PSIP4  is  gebleken  dat  de  resultaten  zodanig  zijn 
dat  de  SIP-meting  en  de  F-PVIS-meting  met  elkaar  vergelijkbaar  zijn.  De  meetaf- 
stand  tussen  het  meetobject  en  bet  explosief  is  echter  wel  groot,  er  zaten  immers 
een  laag  aluminiumfolie,  een  schijfje  PMMA,  een  laagje  siliconenvet  en  twee 
laagjes  siliconenolie  tussen.  Om  nu  te  kijken  of  beide  metingen  met  deze  configu- 
ratie  wel  de  goede  schokgolf-  en  deeltjessnelheid  meten,  is  geprobeerd  dichter  bij 
het  explosief  te  meten. 

Om  dichter  bij  het  explosief  te  meten  is  het  PMMA  bufferplaatje  van  meting  F- 
PVIS4  door  een  laagje  molybdeenfolie  met  een  dikte  van  10  pm  vervangen.  Mo- 
lybdeen  heeft  een  zeer  hoge  dichtheid  zodat  deze  de  schokgolf  nauwelijks  dempt. 
Ook  het  smeltpunt  is  hoog  zodat  eventueel  geioniseerde  lucht  bij  de  springstof  de 
opgedampte  goudlaag  minder  snel  zal  verbranden.  De  rest  van  de  opstelling  blijft 
ongewijzigd. 

4.4.1  De  SIP-meting  (schokgolfsnelheidsmeting) 

Het  resultaat  van  de  SIP-meting  wordt  in  de  twee  onderstaande  tabellen  weergege- 


ven. 
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Tabel  13: 

Resultaten  SIP-meting. 

Lsip  N 

At[s] 

0,978-10-3 

158-10-9 

Tabel  14: 

Berekening 

van  de  druk  in  TNT  uit  de  schokgolfsnelheid  van  de  SIP-meting. 

meting 

Us 

Up 

u  TNT 
Up 

pTNT 

[-] 

[km-s'M 

[krn-s"^] 

[km-s"''] 

[GPa] 

F-PSIP5 

6,190 

2,313 

1,982 

22,979 

4.4.2  De  F-PVIS-meting  (deeltjessnelheidsmeting) 

De  SIP-meting  is  gelukt,  dit  geldt  niet  voor  de  F-PVIS-meting.  Het  niet  verschoven 
lijnenpatroon  is  slechts  aan  een  zijde  van  de  film  goed  te  zien,  aan  de  andere  zijde 
van  de  film  zijn  de  lijnen  vaag  zichtbaar.  Van  het  verschoven  lijnenpatroon  zijn 
maar  aan  een  zijde  drie  punten  te  zien  (zie  bijlage  A).  Na  controle  van  de  F-PVIS- 
opstelling  bleek  dat  deze  was  ingesteld  op  de  interne  reflecties  in  het  systeem.  Dit 
was  de  oorzaak  van  de  het  slechte  beeld  op  de  film.  Daamaast  lijkt  het  dat  de 
goudlaag  minder  reflecterend  is  geworden  tijdens  de  meting.  Er  is  een  deeltjes- 
snelheid  uit  de  film  te  halen  maar  deze  is  waarschijnlijk  niet  goed. 

Bij  deze  meting  is  het  begin  van  de  verschuiving  slechts  aan  een  kant  van  de  film 
te  zien.  Het  niet  verschoven  patroon  is  aan  beide  kanten  redelijk  te  zien  en  kan  dus 
worden  opgemeten.  Met  behulp  van  het  niet  verschoven  lijnenpatroon  kan  het 
midden  van  de  film  bepaald  worden.  Door  nu  de  afstand  tussen  het  verschoven 
patroon  en  het  midden  van  de  film  te  verdubbelen  kan  toch  nog  een  deeltjessnel- 
heid  worden  bepaald.  De  meetfout  wordt  echter  wel  gelijk  tweemaal  zo  groot.  Ook 
het  lichter  en  donkerder  worden  van  de  opname  is  niet  te  zien,  de  lengte  van  de 
film  kan  dus  niet  worden  opgemeten. 
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Tabel  16: 

Berekening  van 
meting. 

de  druk  in  TNT  uit  de  deeltjessnelheid 

van  de  F-PVIS- 

meting 

Up 

u 

Up 

pTNT 

[-] 

[km-s'^l 

[km-s"’] 

[GPa] 

F-PSIP5 

2,862 

2,448 

32,867 

4.5  Gelukte  meting  dichter  op  de  springstof  (F-PSIP6) 


goudlaag 

molybdeenfolie 


glasfiber 


Figuur  /  / ;  Opstelling  bij  meting  F-PVIS6. 

Er  zijn  genoeg  resultaten  van  de  SIP-meting.  Om  de  opstelling  minder  gecompli- 
ceerd  te  maken  wordt  deze  meting  weggelaten.  Er  wordt  geprobeerd  met  de 
F-PVIS  zo  dicht  mogelijk  op  de  springstof  te  meten.  Voor  de  zekerheid  wordt  het 
10  pm  dikke  molybdeenfolie  uit  de  vorige  meting  vervangen  door  molybdeenfolie 
van  30  pm  dik.  Het  aluminiumfolie  deed  alleen  dienst  als  aardelektrode  in  de  SIP- 
meting  en  wordt  dus  weggelaten.  De  siliconeolielaag  tussen  het  goud  en  het  mo¬ 
lybdeenfolie  blijft. 

4.5.1  De  F-PVIS-meting  (deeltjessnelheidsmeting) 

Na  het  ontwikkelen  van  de  film  blijkt  dat  deze  zeer  scherp  en  goed  belicht  is.  Ook 
het  verschoven  lijnenpatroon  is  goed  te  zien  (bijlage  A).  Wat  opvalt  aan  het  ver- 
schoven  lijnenpatroon  is  dat  deze  na  de  stapvormige  verschuiving  nog  een  stukje 
wat  langzamer  oploopt  voordat  deze  zijn  eindwaarde  bereikt.  In  eerste  instantie 
werd  gedacht  dat  dit  de  Von  Neumann-piek  kon  zijn.  Na  het  opmeten  van  de 
lijnafstanden  onder  de  microscoop  blijken  de  verschoven  lijnen  alleen  maar  op  te 
lopen  en  niet  meer  terug  te  lopen,  het  is  dus  geen  piek.  Wat  waarschijnlijker  is  dat 
dit  effect  veroorzaakt  wordt  door  reflecties  in  het  dunne  molybdeenfolie.  Dat  dit  in 
voorgaande  metingen  niet  is  waargenomen  wordt  vermoedelijk  veroorzaakt  door- 
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dat  de  schokgolf  eerst  door  een  andere  materiaal  (molybdeen)  gelopen  heeft  voor- 
dat  deze  het  reflecterend  meetobject  bereikt.  Hierdoor  is  de  deeltjessnelheid  al 
verhoogd.  Als  de  schokgolf  dan  bij  het  grensvlak  met  PMMA  aankomt,  is  er 
slechts  een  kleine  sprong  nodig  om  de  waarde  voor  de  deeltjessnelheid  in  het 
PMMA  aan  te  nemen.  Nu  er  zo  dicht  op  het  explosief  wordt  gemeten  is  het  ver- 
schi!  in  deeltjessnelheid  veel  groter  zodat  deze  meer  tijd  nodig  heeft  om  zijn 
eindwaarde  te  bereiken. 


Op  de  film  is  te  zien  dat  de  lijnen  van  het  verschoven  interferentiepatroon  eerst 
oplopen  en  daama  een  constante  waarde  aannemen.  In  de  volgende  tabel  zijn 
daarom  de  tijdsposities  waarop  de  lijnafstanden  zijn  bepaald  ook  gemeten.  De 
laatste  positie  is  aan  het  einde  van  de  verschuiving,  deze  wordt  gebruikt  voor  de 
tijdberekening  van  de  doorlooptijd  van  de  schokgolf. 

LPMMA.  0,988-10-3  m. 


Tabel  1 7:  Meetresultaten  van  de  F-PVIS  in  meting  F-PVIS6. 


orde 

Di 

[m] 

Dvi 

[m] 

Dv2 

[m] 

Dv3 

[m] 

Dv4 

[m] 

^v5 

[m] 

1 

2,070-10-3 

2,870-10-3 

2,940-10-3 

2,965-10-3 

2,970-10-3 

2,975-10-3 

2 

3,300-10-3 

3,820-10-3 

3,890-10-3 

3,910-10-3 

3,925-10-3 

3,915-10-3 

3 

4,175-10-3 

4,610-10-3 

4,660-10-3 

4,685-10-3 

4,690-10-3 

4,695-10-3 

4 

4,900-10-3 

5,270-10-3 

5,315-10-3 

5,345-10-3 

5,355-10-3 

5,360-10-3 

5 

5,525-10-3 

5,870-10-3 

5,910-10-3 

5,925-10-3 

5,925-10-3 

5,930-10-3 

pos. 

Li 

L2 

L3 

L4 

L5 

Le 

[m] 

[m] 

[m] 

[m] 

[m] 

[m] 

0,000-10-3 

0,100-10-3 

0,200-10-3 

0,300-10-3 

0,400-1 0-3 

1,285-10-3 

Tabel  1 8:  Berekening  van  de  druk  in  TNT  uit  de  deeltjessnelheid  van  de  F-PVIS- 

meting. 


meting 

Up 

u 

Up 

pTNT 

[-] 

[km-s-'l 

[km-s-^] 

[GPa] 

F-PVIS6/1 

1,786 

1,534 

15,075 

F-PVIS6/2 

2,007 

1,722 

18,195 

F-PVIS6/3 

2,107 

1,807 

19,698 

F-PVtS6/4 

2,143 

1,838 

20,253 

F-PVIS6/5 

2,158 

1,850 

20,486 

Tabel  19: 

Berekening  van 

de  druk  in  TNT  uit  de  schokgolfsnelheid 

van  de  F-PVIS- 

meting. 

meting 

Us 

Up 

u  ■'■NT 

Up 

pTNT 

[-1 

[km-s-’l 

[km-s-’' 

[km-s-''] 

[GPa] 

F-PSIP6 

6,151 

2,289 

1,962 

22,578 
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De  posities  uit  label  17  worden  naar  de  tijd  omgerekend  door  deze  te  vermenig- 
vuldigd  met  de  schrijfsnelheid  van  de  streakcamera,  deze  is  125  ns-mm'*.  Om  het 
drukverloop  te  verduidelijken  is  deze  in  een  grafiek  tegen  de  tijd  uitgezet. 


DRUKVERLOOP  IN  F-PVIS6 


tijd  [ns] 


Figuur  12:  Drukverloop  in  F-PVIS6. 


5 


Foutberekening  in  de  F-PVIS  en  SIP-resultaten 


Omdat  de  foutberekening  belangrijk  is  bij  het  meten  en  bij  het  vergelijken  van 
verschillende  meetmethoden  wordt  er  bij  de  foutberekening  uitgegaan  van  de 
resultaten  van  meting  F-PSIP4.  Om  de  drie  meetmethoden  te  kunnen  vergelijken  is 
het  niet  noodzakelijk  dat  de  absolute  fout  telkens  helemaal  wordt  doorberekend  tot 
de  fout  in  het  eindresultaat.  Het  is  voor  een  vergelijking  van  de  methoden  vol- 
doende  om  de  relatieve  fout  van  de  dmk  in  het  PMMA  te  berekenen.  Om  van 
daamit  de  druk  in  het  explosief  te  berekenen  moeten  de  drie  methoden  daama 
dezelfde  berekening  ondergaan.  De  berekende  fout  in  de  druk  in  het  PMMA  die  op 
deze  manier  wordt  gevonden  is  dus  evenredig  met  de  fout  van  de  druk  in  het 
explosief. 

5.1  Fout  in  de  druk  uit  de  met  de  F-PVIS-gemeten  deeltjessnelheid 


De  fout  in  de  deeltjessnelheid,  wordt  berekend  uit  het  veranderende  interferentie- 
patroon  berekend  met  formule  12. 


waarin: 

D,  :  2, 240- 10-3 +  0,0 10- 10-3  m; 

Dv!  :  3, 1 00- 10-3  ±0,01 0-1 0-3  m  Upi  :  2,153-103  ±  53  m-s"'  (2,46%); 

D2  :  3,405-10-3  ±0,010-10-3  m; 

Dv2  :  4,040-10-3  +  0,010-10-3  m  Up2  :  2,216-103  +  73  m-s-'  (3,29%); 

D3  :  4,265  - 1 0-3  ±  0,0 1 0- 1 0-3  m; 

Dv3  :  4,775-10-3  ±0,010-10-3  m  Up3  :  2,161-103  ±  88  m-s-' (4,07%); 

D4  :  5,980-10-3  ±0,010-10-3  m; 

Dv4  :  5,420- 10-3  + 0,010-10-3  m  Up4  ;  2,1 45-103  +  101  m-s-' (4,71%); 

Ds  ;  5, 595 -10-3  ±0,0 10- 10-3  m; 

Dv5  :  5,995-10-3  ±0,010-10-3  m  Ups  :  2,173-103  ±  1 12  m-s-' (5,15%). 

Door  de  kwadratische  verhouding  tussen  Dj  en  Dv  loopt  de  fout  per  orde  op.  Het 
zou  dus,  als  de  kwaliteit  van  de  film  het  toelaat,  het  beste  zijn  om  alleen  de  eerste 
orde  van  de  film  te  gebruiken  bij  de  bepaling  van  de  snelheid.  Om  toch  een  zo 
nauwkeurig  mogelijke  indicatie  te  krijgen  van  de  fout  nemen  we  het  gemiddelde 
van  de  berekende  snelheid  en  fout  over  de  eerste  vijf  ordes,  dit  komt  overeen  met 
het  gewogen  gemiddelde.  De  deeltjessnelheid  is  dan: 

Up  :  2,154  + 0,086  km-s-'  (3,99%). 

De  druk  in  het  PMMA  wordt  berekend  met  formule  13. 
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P  =  (5-Po)-w/4-(Co-Po)-M„  (13) 

De  fout  hierin  is: 

AP  =  (po^™  •  wj)  ■  AS""™  +  (S""™  ■  ul  +  q,""™  •  «p)  •  Ap""™  + 

(pr"  •»,)■  a::,"""  +(2-s-pr“ 

waarin: 

SPMMA  .  1,595+0,015; 

CqPMMA  ;  2,500  ±  0,082  km-s''; 

PqPMMA  .  1 , 1 86  ±  0,050  g-cnr^; 
pPMMA  ;  15,163+ 1,887  GPa(  12,44%). 


5.2  Fout  in  de  druk  uit  de  met  de  F-PVIS-gemeten  schokgolfsnel- 
heid 


Bij  het  berekenen  van  de  schokgolfsnelheid  met  de  F-PVIS  treedt  de  volgende  fout 
op: 


At/.,  = 


^^SIP  ^  ^FILM 


V  ^SIP 


•t/.. 


■'FILM  J 


(15) 


waarin: 

Lhlm  :  1,290-10-3  + 0,005  10-3  m; 
LsiP  :  978-10-6 +  5-10-6  m; 

Us  :  6,065  +  0,054  km-s- '  (0,89%). 


Bij  deze  schokgolfsnelheid  is  ervan  uitgegaan  dat  deze  constant  is  over  de  gehele 
dikte  van  het  SIP-plaatje.  In  praktijk  zal  er  een  demping  optreden.  Hierdoor  is  de 
gemeten  schokgolfsnelheid  lager  dan  de  schokgolfsnelheid  op  het  grensvlak  met 
het  explosief.  Wat  de  bijdrage  is  van  deze  demping  op  de  grootte  van  de  fout  is 
(nog)  niet  in  een  getalswaarde  uit  te  drukken. 


Met  formula  15  kan  de  fout  in  de  deeltjessnelheid,  berekend  uit  de  schokgolfsnel¬ 
heid  met  formula  14,  bepaald  worden: 


t/,,=Q  +  5-M^ 


(16) 


waarin: 

Cq:  plastische  geluidssnelheid  in  het  materiaal  (afhankelijk  van  up)  [km.s‘^  ]; 

S:  materiaalconstante  [-]. 


^A(/ ,  +  AC,  AS  ^ 


t/,.-c, 


■  +  ■ 


(17) 
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waarin: 

Up  :  2,235  ±0,106  km-s'*  (4,74%). 

De  fout  in  de  druk  van  het  PMMA,  berekend  uit  de  schokgolfsnelheid  gemeten  met 
de  F-PVIS,  wordt  berekend  met  formule  12, 

waarin: 

pPMMA  .  16,08  ±2,195  GPa  (13,65%). 


5.3 


Fout  in  de  SlP-resuItaten 


De  fout  in  de  berekening  van  de  schokgolfsnelheid  is: 


AAr 


AC/,  = 

L  M.Jcip  L.\l 

— 

y  ^sip  y 

waarin: 

LsiP  : 

978-10-6  ±5-10-6 

At 

162-10-9  ±2-10-9  s- 

AUs  : 

6,040  ±0,105  km-s-‘  (1,74%). 

(18) 


De  deeitjessnelheid  in  het  PMMA  wordt  berekend  met  formule  1 4,  de  fout  hierin 
wordt  berekend  met  formule  15, 

waarin: 

Up  :  2,219  ±  0,138  km-s-'  (6,22%) 

De  deeitjessnelheid  wordt  omgerekend  naar  de  druk  in  het  PMMA  met  formule  1 1 . 
De  fout  wordt  berekend  met  formule  1 2, 

waarin: 

pPMMA.  15,893  ±2,541  GPa  (15,84%). 
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6  Drukberekeningen  met  de  Kamlet-methode 


De  uitgevoerde  berekeningen  gaan  er  vanuit  dat  de  Hugoniot  van  het  explosief 
bekend  is.  Dit  geldt  echter  niet  voor  alle  explosieven.  Om  deze  berekeningen  toch 
uit  te  kunnen  voeren  wordt  gebruikgemaakt  van  de  Kamlet-methode,  een  empirisch 
bepaalde  benadering  voor  de  drukken  in  explosieve  stoffen.  Voor  TNT  is  die 
weergegeven  in  formule  1 9. 


1,107-po^ +1,702 -PoT  0,655 


(19) 


Eerst  moet  deze  Hugoniot  gespiegeld  worden  om  de  deeltjessnelheid  in  het  explo¬ 
sief  Up™T'.  Samen  met  formule  (16),  waarbij  Us  wordt  vervangen  door  de  detona- 
tiesnelheid,  wordt  deze  omgeschreven  naar; 

P,,  =(2.«™^  -«,)^  •(l,107-p/"’^%l,702.p™"  +0,655)  (20) 

Door  deze  gelijk  te  stellen  aan  de  Hugoniot  van  PMMA  kan  met  behulp  van  een 
spreadsheet  de  deeltjessnelheid  Up™^  berekend  worden.  De  druk  in  de  springstof 
wordt  hieruit  berekend  met  formule  16,  waarbij  Uj  wordt  vervangen  door  de 
detonatiesnelheid.  Hiervoor  moet  dus  op  een  andere  manier  gemeten  worden 
(bijvoorbeeld  met  ionisatiepinnen).  Om  deze  benaderingsmethode  te  kunnen 
vergelijken  met  de  resultaten  uit  het  voorgaande  hoofdstuk  is  tabel  20  opgenomen. 
De  vet  gedrukte  waarde  zijn  de  gemeten  resultaten.  Voor  TNT  is  de  detonatiesnel¬ 
heid  van  6,950  km-s"'  aangenomen. 


Tabel  20:  Resultaten  van  benaderingsmethode  +  samenvatting  resultaten. 


meting 

methode 

Us 

[km-s"’] 

Up 

[km-s"''] 

U  ■'■NT 

Up 

[km-s"'] 

Pberekeninq 

[GPa] 

Pkamlet 

[GPa] 

F-PSIP1 

SIP 

6,300 

2,382 

2,041 

24,129 

22,434 

F-PSIP2 

F-PVIS 

2,640 

2,259 

28,660 

24,656 

F-PSIP2 

F-PVIS 

5,596 

1,941 

1,666 

17,235 

18,615 

F-PSIP2 

SIP 

5,700 

2,006 

1,722 

18,184 

19,180 

F-PS1P4 

F-PVIS 

2,154 

1,847 

20,424 

20,459 

F-PSIP4 

F-PVIS 

6,065 

2,235 

1,916 

21,705 

21,161 

F-PSIP4 

SIP 

6,040 

2,219 

1,903 

21,455 

21 ,025 

F-PSIP5 

F-PVIS 

2,862 

2,448 

32,867 

26,595 

F-PSIP5 

SIP 

6,190 

2,313 

1,982 

22,979 

21 ,838 

F-PSIP6 

F-PVIS 

1,786 

1,534 

1 5,075 

17,268 

F-PSIP6 

F-PVIS 

2,007 

1,722 

18,195 

19,186 

F-PSIP6 

F-PVIS 

2,107 

1,807 

19,698 

20,052 

F-PSIP6 

F-PVIS 

2,143 

1,838 

20,253 

20,364 

F-PSIP6 

F-PVIS 

2,258 

2,158 

20,486 

21 ,359 

F-PS1P6 

F-PVIS 

6,199 

2,319 

1,987 

23,072 

21,886 
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Conclusie  en  discussie 


Uit  meting  F-PSIP4  blijkt  dat  de  resultaten  van  de  drukmetingen  volgens  de  SIP- 
methode  en  de  F-PVIS-methode  goed  met  elkaar  overeenkomen. 

Toch  is  drukbepaling  uit  de  deeltjessnelheid  met  de  F-PVIS  de  beste  methode.  Bij 
een  goede  film  kunnen  de  lijnafstanden  met  een  fout  van  0,01  mm  opgemeten 
worden.  De  fout  in  de  berekening  van  de  uiteindelijke  druk  is  dan  ongeveer  12,5% 
en  wordt  voomamelijk  veroorzaakt  door  de  diverse  omrekeningen.  Bij  een  on- 
scherpe  film  worden  de  lijnen  dikker  wat  het  eindresultaat  sterk  beinvloed,  een 
goede  film  is  dus  een  voorwaarde. 

Als  de  druk  volgens  de  SIP-methode  wordt  bepaald  is  de  fout  hierin  bijna  16%. 
Deze  grote  fout  wordt  vooral  veroorzaakt  door  de  fout  die  optreedt  bij  het  aflezen 
van  de  tijd  van  de  schokgolfdoorgang  op  de  DSA.  De  relatief  grote  stijgtijd  van  het 
signaal  bij  het  intreden  van  de  schokgolf  in  het  plaatje  zorgt  voor  een  meetfout  van 
5  ns.  Verder  moet  de  hieruit  afgeleide  schokgolfsnelheid  eerst  omgerekend  worden 
naar  de  deeltjessnelheid,  dit  veroorzaakt  een  extra  fout.  Een  voordeel  is  dat  het 
signaal  wat  gemeten  wordt  ook  echt  het  in-  en  uittreden  van  de  schokgolf  in  het 
SIP-plaatje  is,  echt  extreem  grote  afwijkingen  kunnen  dus  niet  optreden. 

Dit  geldt  niet  als  de  druk  uit  schokgolfsnelheid,  gemeten  met  de  F-PVIS,  berekend 
wordt.  Het  is  ook  een  goede  methode  maar  er  is  een  onzekere  factor.  Het  is  name- 
lijk  moeilijk  te  bepalen  of  de  grinlens  goed  tegen  het  PMMA-plaatje  aanzit.  De 
dikte  van  het  plaatje  wordt  met  een  fout  van  ±0,005  mm  opgemeten  en  als  de 
grinlens  niet  goed  tegen  het  PMMA-plaatje  aanzit  wordt  deze  fout  snel  groter.  Dit 
heeft  een  grote  invloed  op  de  fout  in  de  berekening  van  de  schokgolfsnelheid  en 
dus  ook  op  de  berekening  van  de  uiteindelijke  druk  in  het  explosief.  Bij  een  goede 
opstelling  is  de  fout  in  de  berekende  druk  uit  de  schokgolfsnelheid  bepaald  met  de 
F-PVIS  ongeveer  1 3,5%.  Dit  is  beter  dan  de  SIP-meting  maar  aan  het  resultaat 
moet  dus  in  vergelijking  met  de  SIP-meting  minder  waarde  worden  gehecht.  Een 
groot  voordeel  van  deze  meetmethode  is  wel  dat  de  gegevens  voor  berekenen  van 
de  schokgolfsnelheid  altijd  bij  de  gegevens  van  de  deeltjessnelheidsmeting  zitten. 
Met  deze  methode  kunnen  dus,  mits  de  opstelling  goed  is,  twee  schokgolfgroothe- 
den  in  een  meting  bepaald  worden.  Hiermee  is  een  punt  uit  de  Hugoniot-curve  van 
het  explosief  bepaald,  normaal  zou  hiervoor  volgens  twee  methoden  gemeten 
moeten  worden. 

De  benaderingsmethode  voor  de  Hugoniot  van  een  springstof  werkt  goed.  Het 
verschil  tussen  de  berekende  en  benaderde  waarde  van  de  druk  ligt  binnen  de  fout 
van  deze  druk. 

Het  PMMA  heeft  een  dempende  werking  op  de  schokgolf.  Dit  is  te  zien  aan  de 
afname  van  de  deeltjessnelheid,  en  dus  van  de  schokgolfsnelheid,  in  de  film  van 
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meting  F-PSIP2.  Door  deze  dempende  werking  kan  er  worden  geconcludeerd  dat 
voor  een  schokgolfsnelheidsmeting  het  PMMA-plaatje  niet  te  dik  mag  zijn,  er 
wordt  immers  een  gemiddelde  snelheid  bepaald. 

Het  plaatje  in  meting  F-PSIP2  van  3,94  mm  is  duidelijk  te  dik,  beide  schokgolfme- 
tingen  vallen  met  respectievelijk  5,70  km-s‘*  en  5,59  km-s'*  uit  de  toon  met  andere 
gemeten  schokgolfsnelheden  die  alien  boven  de  6  km-s"'  liggen.  De  gemeten 
schokgolfsnelheden  in  meting  F-SIP4,  waar  een  1  mm  dik  PMMA-buffer  is  ge- 
bruikt  en  meting  F-PVIS6,  waar  alleen  een  stukje  molybdeenfolie  van  30  pm  dik 
op  het  TNT  is  geplaatst,  verschillen  niet  veel.  De  1  mm  dikke  PMMA-buffer  heeft 
dus  geen  grote  invloed  op  het  resultaat.  Een  grote  afwijking  komt  dus  voor  als  er 
meer  dan  2  mm  PMMA  gebruikt  wordt.  Desondanks  wordt  er  toch  naar  gestreefd 
zo  dicht  mogelijk  op  het  explosief  te  meten,  dit  is  altijd  beter. 

De  meest  ideale  opstelling  om  de  druk  in  een  explosief  te  meten  is  de  opstelling 
van  meting  F-PVIS6  en  maakt  gebruik  van  30  pm  dik  molybdeenfolie.  Pogingen 
om  nog  dichterbij  het  explosief  te  meten  met  10  pm  dik  aluminiumfolie  (F-PVIS2) 
en  10  pm  dik  molybdeenfolie  mislukten  door  een  slecht  zichtbaar  verschoven 
lijnenpatroon  op  de  film.  Ook  een  herhaling  van  meting  F-PVIS5  met  een  nieuwe 
elektronische  streakcamera  (gegevens  hiervan  zijn  helaas  nog  niet  beschikbaar)  liet 
een  vaag,  niet  meetbaar,  strepenpatroon  zien.  Mogelijk  is  op  deze  kleine  afstand, 
de  hitte  veroorzaakt  door  het  explosief  zo  intens  dat  de  goudlaag  wordt  aangetast. 
Een  andere  mogelijkheid  is  dat  de  grote  druk  die  in  de  Von  Neumann-piek  heerst 
toch  iets  buiten  het  explosief  doorwerkt  en  de  goudlaag  aantast. 
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8  Verklaring  van  gebruikte  symbolen  en  woorden 


AX :  Golflengte  verschuiving  door  Doppler-effect  [m] 

At :  Doorlooptijd  van  de  schokgolf  door  het  SIP-plaatje  [s] 

X :  Golflengte  licht  [m] 

(p :  Faseverschil  [-] 

Po’^,  Pi’^ :  Soortelijke  massa  voor  en  na  de  schokgolf  in  medium  x  [kg-m'^]* 
0 ;  Hoek  tussen  invallende  lichtstraal  en  normaal  [rad] 

c :  Lichtsnelheid  [m-s  *] 

Co :  Plastische  geluidssnelheid  [m-s"'] 

d  ;  Afstand  tussen  spiegels  van  Fabry-Perot  etalon  [m] 

D :  Detonatiesnelheid  [m-s‘*] 

Dj :  Diameter  van  ring  j  [m] 

Dj+i  :  Diameter  van  ring  j+ 1  [m] 

Dv ;  Diameter  van  verschoven  ring  j  [m] 

Eo’^,  E|’‘:  Inwendige  energie  van  medium  X  voor  en  na  de  schokgolf  [J-m’^J 

Fc :  Fringe-constante  [m-s"*] 

:  Hugoniot  van  medium  x  [-] 

Ij ;  Lichtintensiteit  van  doorgelaten  licht  [W] 

;  Lichtintensiteit  van  gereflecteerde  licht  [W] 

fj. ;  Brandpuntsafstand  van  lens  C  [m] 

Lfiim  :  Lengte  van  de  verschoven  lijnen  op  de  film  [m] 

Lsip :  Dikte  van  het  SIP-plaatje  [m] 

P :  Druk  [Pa] 

Pq  :  Chapman-Jouguet-druk  [Pa] 

Pq^,  Pl’^:  Druk  voor  en  na  de  schokgolf  X  [Pa] 

R ;  Reflectiecoefficient  [-] 

S  :  Constante  [■] 

t9o% :  Fill-time  [s] 

Up*:  Deeltjessnelheid  in  medium  X  [m-s"'] 

Up:  Deeltjessnelheid  [m-s"’] 

Us:  Schokgolfsnelheid  [m-s''] 

v:  Snelheid  [m-s-'] 

Vq’^,  Vi’^:  Inverse  soortelijke  massa  [m^-kg"'] 

Z:  Schokgolfimpedantie  [kg-m^-s"'] 


In  formules  die  met  een  ‘curve-fit’  zijn  gemaakt,  hebben  de  gemeten  groothe- 
den  de  eenheid  [GPa],  [km-s"']  of  [kg-m"^]. 
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Figuur  A.l :  Film  van  meting  F-PSIP2,  reflectie  met  een  goudlaag. 


Figuur  A.l :  Film  van  meting  F-PSIP4,  schokgolfdemping  met  PMMA. 


J, 


Figuur  A.3:  Film  van  meting  F-PSIP6,  dichter  op  de  springstof. 
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